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Abstract 

We introduce the concept of natural Poisson bivectors, which generahzes the Benenti 
approach to construction of natural integrable systems on Riemannian manifolds and 
allows us to consider almost the whole known zoo of integrable systems in framework of 
bi-hamiltonian geometry. 



1 Introduction 

Let us consider a natural integrable by Liouville system on a Riemannian manifold Q for 
which the Hamilton function 

n 

H = T + V=Y, &jmPj + Viqi,...,qn) (1.1) 

is the sum of the geodesic Hamiltonian T and potential energy V. Integrability means that 
there are functionally independent integrals of motion Hi = H, H2 , ■ ■ ■ , Hn in involution 

{Hi,Hj} = {dHi,PdHj) = 0, 

with respect to the canonical Poisson brackets defined by the following Poisson bivector 



Among integrable systems we want to pick out a family of bi-integrable systems for which: 

• There is a second Poisson bivector P' compatible with P, i.e. 

[P,P'] = 0, [P',P'] = 0, (1.3) 
where [., .] is the Schouten bracket. 

• Integrals of motion Hi, . . . , Hn are in bi-involution with respect to both Poisson brackets 

{H„Hj} = {H„Hjy = 0. (1.4) 
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First example of bi-integrable systems are bi-hamiltonian systems. The concept of bi- 
hamiltonian vector fields was introduced firstly by Magri studying tlie Korteweg-de-Vries 
equation in order to explain integrability of soliton equations from the standpoint of clas- 
sical analytical mechanics [28]. However, for the overwhelming majority of known natural 
integrable systems on Riemannian manifolds the Hamiltonians H (jl.ip give rise to non bi- 
hamiltonian vector fields X = PdH. The natural obstacle for existence of the bi-hamiltonian 
vector fields in finite-dimensional case is discussed in [lOj. 

The second special but more fundamental example of bi-integrable systems are separable 
systems, for which there exist n separation relations of the form 



(j)i{ui,Vi,Hi,...,Hn) = , i = l,...,n, with det 



/O. (1.5) 



Here u = (ui, . . . , u„) and v = {vi, . . . , Vn) are canonical variables of separation, {uj, Vj} = 6ij. 
The proof of the fact that any separable system is a bi-integrable system with respect to second 
Poisson bivector 

P' = ( n ) ' F = diagf/i(ni,ui), . . . ,/„(n„,t;„) 



^ -F ^ 

labeled by n arbitrary functions fk on variables of separation may be found in [39]. The main 
problem is how to describe these bivectors in terms of initial physical variables or, equivalently, 
how to determine variables of separation for a given natural Hamilton function on Q. 

This problem is solved for the quadratic in momenta integrals of motion Hi,...,Hn 
separable in orthogonal coordinate systems on Q, see O [3l HI [U EJ [HI [El [2ll [22l [35] and 
references within. In this case initial physical coordinates are related with coordinates of 
separation by the point transformation 

n 

Qi = Qiiui, . . . ,Un), Pi = '^hij{ui, . . . ,Un)Vj . (1.6) 

i=i 

The first algebraic condition for separability systems with quadratic in momenta integrals of 
motion has been found by Stackel [35] . Then in [26] Levi-Civita proved that a Hamilton- Jacobi 
equation H(q,p) = E admits a separated solution if and only if the separability conditions or 
separability equations of Levi-Civita are identically satisfied 

diHdjHd'd^H + d'd^HdidjH - diHd^Hd'djH - d'HdjHdid^H = , i,j = l..n, 

here dk = d/duk and = d/dvk ■ Using transformations ()1.6p it is easy to rewrite Levi- 
Civita equations as polynomial equations of fourth degree in the momenta pi, . . . ,p„ and to 
note that fourth-degree homogeneous part of the Levi-Civita equations depends only on the 
geodesic Hamiltonian T. Such as Levi-Civita equations must be identically satisfied for all 
admissible values of pi, . . . ,p„ it means that: 

• the separation of the geodesic equation is a necessary condition for the separation at 

• the study of the geodesic separation plays a prominent role. 

As was shown by Eisenhart (for an orthogonal case at gjj = for i ^ j) and by Kalnins and 
Miller for generic case, the geodesic separation is related to the existence of Killing vectors 
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and Killing tensors of order two [23} [23] . These ideas are nicely embraced by the geometric 
notion of Killing webs discussed in [21 [3l [3] • 

In framework of the geometric Benenti theory it is possible to construct a basis of Killing 
tensors by coordinate independent algebraic procedure starting with special tensor L with the 
following properties: 

1. L is a conformal Killing tensor of gradient type, 

2. the Nijenhuis torsion of L vanishes, 

3. L has pointwise simple eigenvalues. 

According to [6l [IHl [22l |37] , these conditions entail that we can define the second Poisson 
bivector 



I 



\ 







-L 



E 



dLki dL 



\ 



kj 



dqj dqi 



Pk 



(1.7) 



and that eigenvalues of the recursion operator = P'P~^ are the desired variables of sep- 
aration. Some algorithms and software for calculation of the Benenti tensor L starting with 
a given natural Hamilton function on the Riemannian manifold Q of constant curvature may 
be found in [181EI1I49]. 

In this note we consider natural Hamiltonians (jl.ip on M^" with unit metric tensor 



H = '^Pi +V{qi,...qn) , gi. 



1=1 



and the corresponding natural Poisson bivectors 



P' 



k=l 



dpi dpj 



Qk Ilj 



V 



+ 



/ 



V 



A, 



-Aji 



E 

k=l 



dA 



ki 



dAkj 

dqj dqi 



Pk 



(1.9) 



which are separable on geodesic and potential parts. Here geodesic Hamiltonian T and geodesic 
matrix H depend only on momenta pi, . . . ,pn and H has zero Nijenhuis torsion as a tensor field 
on M" with these coordinates. Potential V and potential matrix A depend only on coordinates 
qi, . . . ,qn and the Nijenhuis torsion of A on configurational space Q = M" is equal to zero. 

This paper belongs mostly to so-called "experimental" mathematical physics. Our main 
aim is to present some of the most striking examples of matrices H, A and to discuss the yet 
unsolved problems: 

• invariant definition and construction of H and A without an unpretentious solution of 
the equations (jl.3p in the special coordinate system fixed by unit metric tensor (jl.Sp : 



finding the integrals of motion in bi- involution from P'; 
calculation of separation variables; 

definition of similar bivectors on another Riemannian manifolds. 
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Note that at 11 = we have exactly the same problems, which have been partially solved 
in framework of the Benenti theory. In fact the Benenti recursion procedure concerns with 
a special subclass of the Stackel systems, which contains for instance the classical Jacobi 
separation of the geodesic flow on an asymmetric ellipsoid [U [5] . Another solutions of these 
problems are related with lifting of the Stackel systems to bi-Hamiltonian systems of Gelfand- 
Zakharevich type using extension of the initial phase space [6l \T9\ I22j . However, all these 
solutions are relatively simple because momenta pj are linear functions on vj (I1.6P and, there- 
fore, our Hamilton function H has natural form simultaneously in physical variables {p, q) 
and in variables of separation {u,v). 

In generic case at H 7^ for the present we do not have any satisfactory solutions of these 
problems. Some efforts of applying the Killing theory to integrable systems with higher order 
integrals of motion have been made in |30] . Particular solutions of the other problems by 
brute force method may be found in [27113811111 H2l H5l[i6lli8] . 

The paper is organized as follows. In Section 2 the concept of natural Poisson bivectors 
on the Riemannian manifolds is briefly reviewed. The Toda lattices and rational Calogero- 
Moser systems illustrate possible applications of this concept. In Section 3 the problem 
of classification bi-integrable systems on low-dimensional Euclidean spaces is treated. The 
Henon-Heiles system, the system with quartic potential, the Holt-like system and some new 
bi-integrable systems are discussed as well. In Section 4 we introduce natural Poisson bivectors 
on the sphere S". At n = 2 we show natural bivectors associated with the Kowalevski top, 
Chaplygin system and Goryachev-Chaplygin top. The possible generalizations of natural 
Poisson bivectors are discussed in final Section. 



2 Natural Poisson bivectors on the Riemannian manifolds 



Let Q be an n-dimensional Riemannian manifold. Its cotangent bundle T*Q is naturally 
endowed with canonical invertible Poisson bivector P, which has the standard form (|1.2p in 
fibered coordinates {p,q) on T*Q. 

Definition 1 A Poisson bivector P' on T*Q has a natural form if it is a sum of the geodesic 
Poisson bivector P^ and the potential Poisson bivector defining by torsionless (1,1) tensor 
field A{qi, . . . ,qn) on Q 



( 



P' = Pt + 



\ 



-A 



E 

k=l 



dA 



ki 



dA 



kj 



dqj dqi 



Pk 



(2.1) 



In fact, here we assume that bi-integrability of the geodesic motion is a necessary condition 
for bi-integrability in generic case at V ^ and, therefore, natural bivector P' has to remain 
the Poisson bivector at A = 0. It allows us to classify all possible geodesic solutions P^ of the 
equations (II. 3p and then to add to them various consistent potential parts. 

The main result of the our experiments is that all the known geodesic Poisson bivectors 
P^ have the following form 

dpU H 

Pt={ I • (2.2) 

-H^ dgU 
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Here entries of geodesic matrix 11 are the homogeneous second order polynomials in momenta 



n 



(2.3) 



k,m=l 



similar to the geodesic Hamiltonian T (jl.lh , whereas antisymmetric tensors 9^11 and 9^11 are 
given by (jl.9p or by more complicated expressions considered in Section 12.2^ Section [5] and 
Section 15. 1[ 

Proposition 1 Natural Poisson bivectors P' (KJSKM '^'^^ unambiguously determined by a 
pair of geodesic and potential n x n matrices (H, A) on 2n- dimensional space T*Q. 

Our definition of the natural Poisson bivectors P' drastically depends on a choice of coordinate 
system. We hope that further inquiry of invariant geometric relations between metric tensor 
G and matrices 11, A allows us to get more invariant and rigorous mathematical description 
of these objects. 

There is a family of integrable potentials V, which may be added to a given geodesic 
Hamiltonian T in order to get integrable Hamiltonian H = T + V. By analogy, we define a 
family of compatible potential matrices A, which may be added to a given geodesic matrix H 
in order to get natural Poisson bivector P' compatible with P. 

Definition 2 Potential matrix A is compatible with geodesic matrix H if the natural Poisson 
bivector P' ^2.1'^2.S^ satisfies the equations I so that P — XP' is a Poisson bivector for 
each A. 

Remark 1 In some sense we have to waive previous the principle that geodesic motion is 
completely independent from potential. In fact for the same geodesic Hamiltonian T we have 
many different geodesic matrices H compatible with various potential matrices A, i.e. our 
choice of H depends on potential V. 

2.1 Integrals of motion versus variables of separation 

Below we suppose that natural Poisson bivector P' (12. ip is always compatible with canonical 
bivector P, so that the phase space T*Q becomes bi-hamiltonian manifold endowed with the 
hereditary recursion operator 

N = P'P-^ . 
In general, there are three different occasions: 

1. recursion operator produces the necessary number of integrals of motion; 

2. recursion operator generates variables of separation instead of integrals of motion; 

3. recursion operator produces only some part of the integrals of motion or variables of 



In the first case, the traces of powers of the recursion operator are functionally inde- 
pendent constants of motion 



separation. 




Hk = —trace N'' , k = 2, . . . 



n. 



(2.4) 
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and our natural Hamiltonian H (jl.ip is directly determined by 11 and A: 



H = Hi, T = ^nii, and V = Y,^ii- 

i=l i=l 

This Hamiltonian defines natural bi-Hamiltonian system on T*Q. 

In the second case, the traces of powers of the recursion operator N remain functionally 
independent constants of motion for an auxiliary bi-Hamiltonian system on T*Q, which differs 
for our target system with Hamiltonian H (II. ip 

^ n n 

Hy^-traceN, Ty^^Ua, and V^J^^a- 

i=l i=l 

In this situation we will treat eigenvalues Uj of 

B(A) = (^det(iV- AI) j = (A-ni)(A-n2)---(A-'u„), (2.5) 

as separation variables for a huge family of separable bi-integrable systems on T*Q associated 
with various separated relations (II. 5p . Of course, this construction will be justified only if we 
are capable to obtain desired Hamilton functions H p.ip from (jl.5p . 

In third case recursion operators N is degenerate and produces only a part of integrals of 
motion or variables of separation. We can leave this difficulty and get the necessary number 
of functionally independent integrals of motion using some additional assumptions about the 
so-called control matrices F defined by 

n 

P'dH, = P^FijdHj, i = l,...,n. (2.6) 

According to [38], we can fix some special forms of F and try to get the corresponding natural 
Hamiltonians Hi and natural Poisson bivectors P' simultaneously. 

Now we are going to illustrate the first and the third opportunities by examples of the 
n-body Toda lattice and of the rational Calogero-Moser system, respectively. Construction of 
separation variables are considered in Section 13.31 and Section [5l 

2.2 The Toda lattices 

Let us start with a well-known second Poisson tensor for the open Toda lattice associated 
with An root system 



Til r\ r\ Tl' r\ r\ r\ r\ 

P Qi-g^+i A A A f2 7) 

dPi+i dpi ^^'dqi dpi f^- dqj dqi' 

This bivector is different from the natural Poisson bivector defined by (j2.m2.2p . Nevertheless, 
using recursion operator = PP~^, it is easy to get a quadratic in momenta Poisson bivector 

P' = NP, (2.8) 
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which has the natural form of (j2.m2.2p if we put 



U = diag{pl,...,pl) 



dgU = 0, 



1 fdUii dUjj 



+ 



2 V dpi dpj 
0, 

1 /an,v du 



J3 



i = j; (2.9) 

i > j; 



2 V dpi dpj 

and if n X n potential tensor A = —E A is a product of two antisymmetric matrices with 
entries 

1, i < j; 

0, i = j; and = a e«"-'?'+i . (2.10) 

-1, i>j; 

In this case the recursion operator = P'P~^ produces the necessary number of integrals of 
motion Hk defined by (j2.4p and the Hamilton function has the natural form 

n 71—1 

Hi = - trace = ^ + 2a ^ ^'^-Q^+i . 

i=l i=l 

In order to get variables of separation we have to introduce another linear in momenta Poisson 
bivector [39], which may be rewritten in natural form as well. 

Remark 2 In fact matrix A ()2.10p is the well-known second matrix in the Lax equation 
L = [L, A] for the open Toda lattice. In framework of the group theoretical settings of 
integrable systems the Lax matrices are viewed as a coadjoint orbits of a Lie algebras. We 
believe that natural bivector P' (j2.8p has transparent algebro-geometric justification, similar 
to compatible bivectors from [7j. 

The Poisson bivectors for the Toda lattices associated with BCn and P„ root systems 
have the natural form of (12.1112.21) if 



n = diag(p2,...,p2)^ dgU = 0. 



dpi 
0, 



dpi 



i = i; 



and if n X n potential parts are given by 

BCn A = -{l + 2E)A + be'i" B + ce'^'i" C 



(2.111 



A = -{1 + 2E)A + dei"-^+i" D, b,c,d,G 



where I is the unit matrix, A is given by (j2.10p and E is a strictly upper triangular matrix 



Ei, 



1 i<j 
i>j 



Matrices B, C and D have non-zero entries only in the last columns: 

/ 



B 



( 


V 



1\ 

1 
1/ 



c 





Vo 



2 \ 

2 
1/ 



D 






vo 



2 2 \ 

2 2 

2 1 
0-10/ 
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As above, recursion operators N generate integrals of motion Hi^. (j2.4p and Hamilton functions 
have the natural form 

n n—1 

2(Jn 



ce 

i=l i=l 



n-1 



i=l i=l 



These natural Poisson bivectors in (p, g)-variables have been obtained in [12], whereas bi- 
hamiltonian structures for the periodic Toda lattices and construction of the separation vari- 
ables are discussed in [lOl [l3] . 



Remark 3 The " relativistic" modification of the natural Poisson bivector (j2.ip associated 
with relativistic n-body Toda model [31j is considered in Section (j5.2p . 

2.3 The Calogero-Moser system 

The bi-Hamiltonian formulation of the Calogero-Moser system can be found in [H [T6l I29j . 
We present new and very simple natural Poisson bivector (jl.9p . which is different form these 
known Poisson brackets expressed directly in terms of integrals Hk ■ 

The n-particle rational Calogero-Moser model associated with the root system An is 
defined by the Hamilton function 

1=1 

where a is a coupling constant. The second natural Poisson bivector P' (jl.9p for this system 
is defined by symmetric geodesic matrix 

Yl=p®p, Uij=piPj (2.13) 

and potential matrix A with entries 

n 2 

In this case recursion operator N = P'P^^ generates only the Hamilton function 

traceA^'' = 2(2i?)'' and PdH = P'dlnH. (2.15) 



and we have so-called irregular bi-Hamiltonian manifold. Nevertheless, integrals of motion 
1 

Ik 



Hk = — trace obtained from the standard Lax matrix 



kl 

"1 01—09 Ot—O^ 

£ _ (J2-IJ1 



91-92 11- Qn 

P2 



... p„ 
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are in bi-involution (jl.4p with respect to the Poisson brackets defined by P (|1.2p and P' (j2.13t - 
I2.14p . Besides these n integrals of motion the rational Calogero-Moser system admits n — 1 
additional functionally independent integrals of motion 

Km = mgiHm - gmHi, gm = ^ ^^^^Qj^ , Hm^ , 171 = 2, ...,71. 

All these integrals of motion Hi^. and Km may be obtained from the Hamilton function H = H2 
(|2.12p as polynomial solutions of the following equations 

PdH=^P'dlnHk = ^—P'dlnKm. (2.16) 
k m — 1 

In Section 13.21 we discuss solutions of the equations (|2.16p associated with another natural 
Poisson bivectors and other bi-integrable systems. 

Remark 4 We suppose that trigonometric (elliptic) Calogero-Moser systems and their gener- 
alizations associated with other root systems may be associated with natural Poison bivectors, 
see [38] at n=2 . In order to describe the Ruijsenaars-Schneider model we can try to introduce 
another Poisson bivector P' similar to relativistic Toda case, see Section 15.21 



3 Natural bivectors on low-dimensional Euclidean spaces. 

Now let us come back to the Poisson manifold i?^" endowed with the natural Poisson bivector 
(jl.9p . In this case the corresponding Poisson bracket {.,.}' looks like 



9fc, 



{qi,Pj} = Uij + Aij , {qi, qj} = ^ - j 

The geodesic Hamiltonian T is second order homogeneous polynomial in momenta, so we are 
going to suppose that entries of 11 are second order homogeneous polynomials in momenta as 
well. This assumption allows as to get a lot of natural Poisson bivectors P' (II. 9p compatible 
with canonical bivector P and describe the corresponding bi-integrable systems. For brevity 
we will not consider the complete classification and restrict ourselves by discussing only more 
interesting examples for n = 2, 3. Some other examples may be found in \n\ I38j. 

3.1 Integrals of motion via recursion operator 

In this section we suppose that geodesic Hamiltonian T and potential V are directly deter- 
mined by H and A (p^ . 



Case 1: Let us start with the non-degenerate geodesic matrix 



= 1 ( + i^"^ ° ). BO that r=?i±M. (3.2) 
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There are some potential matrices A compatible with it, for instance: 

: — Oi + Cl On + C9C/0 + C^Qo — H : 



Ad) 



ciql 



32 



16 +^^ + ^2g2 + |)gi^ 



Ciql + 02^2 + C3g2 



^^'^ = |- 92(3^? + 16gi) + C2 (2^2 + f ) + 20392 



and 



A(i') 



^ + (3ci9^ + 02) — + ciq^ + C2qi + ^ — ^ — + \ 2ciq^ + 02^2 



(3.3) 



C3\ \ 



ciqfq2 



4 I 2 I ''3 
92 



= f (gt + 6,?,i + 8^2^) + + 4gi) + ^ . 

4 ^52 



(3.4) 



The second integrals of motion H2 = trace N are fourth order polynomials in momenta: 

.2/Q^,„_ I ^ „3 ^4 

il 

2 2 "'"^ 32 

and 

-,4 



[r(l) 4 , 9l(3cl92 + C2) 2 Cigf 9l/ 2/« 2 , 2x , /n ,4 \ O 2\ /o r-^ 

=Pi + Pi —piP2 -^[ci[Qq2 +qi) + ci(8c3 + 4c2g2) - 2c2 j (3.5) 



H. 



(1') 



(pI + C2qlf + CiqlUql + QqDpl + qfpj " Mmpm) + ^^^^ (3.6) 

^ ^ qi 

+ \ciqj{qj + 2ql){ciql + 2ciql + 4c2) • 

These integrable systems were found by using the weak-Painleve property of equation of 
motion and the direct search of fourth order polynomial integrals of motion, see jl7| I20j. 

Remark 5 The Henon-Heiles system with potential V^^'^ and the system with fourth order 
potential V^^"> admit various integrable generalizations |17^ I20j. which are considered in the 
Section 



Case Ig: At n = 3 the immediate generalization of H^^) (|3.2p looks like 



n(i9) 








3 2 




PlP2 


pI 


2 


2 


V2piP2 


V2pl 


V 2 


3 



V2piP2 2piP3 \ 

3 3 

V2pl 
6 

¥-3 ; 



P1+P2+ pI 



(3.7) 
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One of the potential matrices A compatible with H^^^) is equal to 



Ads) 



ci 



q2{9ql+?,ql) 
2 



99? 



4 



4g: 



8 

3V2(7f(g2+v^g3) 



3V^..(3,^+4,H2,g+272,293) 372,^(4,2+72,3) 3^2^^ + 3^^2^3 + ^^3 ^ 



+ C2 



2 yi 



+ C3 



V 6(^/2(72 + q^)qi ^V^ql ^ql + ^q\ j 



/ „ 91 \/2gi 

g2 

3\/2gi |-i \/2g3 



SO that 



ci 



?i + ^qlq2 + sV2q3ql + V2qi] + C2(16gi + + 4ql) + C3 f 2^2 + ^ 



93 



In this case the integrals of motion H2 = trace N"^ and = det are the fourth and sixth 
order polynomials in momenta, such that few pages are necessary to write them down. 



Remark 6 At 02,3 = this potential is equivalent to potential Vio in [32]. Similar Poisson 
bivectors may be constructed for another potentials from |32J and for n dimensional general- 
izations of the Henon-Heiles systems and systems with quartic potentials |14j . 

3.2 Integrals of motion via control matrices 

In [38] we have obtained some natural Poisson bivectors using the Lenard and Probenius 
control matrices. Now we present some new examples of two-dimensional bi-integrable systems 
associated with degenerate control matrix 



Hi 



(3.8) 



It means that Hamiltonian Hi is the solution of the equation 

X = PidHi = P'd In Hi, 
whereas H2 is the solution of the equation depending on rational parameter x 

X = PidHi = xP'dlnH2, 



(3.9) 



similar to the equations for the Calogero-Moser system (|2.16p . 

An algebraic construction of such two-dimensional Hamilton functions Hi, additional 
integrals of motion H2 and natural Poisson bivectors P' (jl.9p has been proposed in [27]. 

Case 2: Let us consider the degenerate symmetric matrix associated with the Calogero-Moser 
systems (I2.13j) 



n(2) = 1 ( Pi P^P^ 
2 V P1P2 pI 



such that 



T 



P1+P2 



(3.10) 
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There are some other potential matrices A compatible with it, for instance, 

■v2 



A(2) 



Cl 



H 



(2) 



(9? + ql? 

T + , 



91 9192 

9192 92 
r(2) 



y(2) 



Cl _ Cl 

9i + 92 '^^ 



(3.11) 



and 



A(2') 



/ Cl (d + 2)c29f 
9? 2gf 2 



dC292"^ \ 



ci92 _^ + 2)c2g2 



d+l 



H 



(2') 



V 9? 

T + y(^) 



9f+^ 



2gf-3 
dc2qi 



1/(2') ^ £l , ^^2 92 

9? ^ 9^^ 



(3.12) 



H. 



(2') 



bi92-P29i)^ + ^ + (-l)^-^^^^^%i^') <^ 

9i 9i y 



In both cases second integrals of motion H2 were found as the solutions of the equation (|3.9 
at X = 1. Let us note that the matrix A^^'^ (I3.12p is a particular case of matrices fixed by 



1 



92 



A12 = ^ $ — 
9? V91, 

where ^ is an arbitrary function [38]. 

Case 2g: System (j3.10l - l3?TT]) has an obvious n-dimensional counterpart 



n 



(2s) _ PiPj 



(2g) _ cig^gj 



2\2 ' 



Cl 



i=l 



(3.13) 



Remark 7 According to ^J, at n = 2 there are a lot of integrable deformations of the 

Cl 

centrally symmetric potential ^ = for which the forth and the sixth order polynomial 

integrals H2 do not the products of the Hamilton function H and the second order polynomials 
as in (13. lip . It will be interesting to study similar deformations at n > 2. 

Case 3: Now let us consider another metric and degenerate non-symmetric matrix 

"'''=Kt;|^ .pi)- r=(«+%ip.. (3.14) 

It is easy to prove that matrix A may be added to II^^^ if and only if 

_2b _b 

A12 = 91 $(929i 
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For instance, if ^{z) = z'^ and 7 = (26 + a + bd) then 



A{3) 



7 



a + b 



\ a[a + 0) a 



aci{d+l) . -M£*±2) 



a + 6 



■qi qi 



bci{d + l) a+i - 
' a + b «i 



T^f3) d+1 I -(a+6)/a 

y^'''' = Ci " ^2 + C2 9i 



(3.15) 



In order to get an additional integral of motion H2 we have to use the equation p.Op because 
the recursion operator is degenerate. Depending on values of a, b and d second integral 
of motion H2 may be second, forth or sixth order polynomial in momenta. For instance, we 
present some examples with forth order polynomial integrals of motion 



H2 = Ap\ - 10 {?>piql - 30^1^29192 + 2bplql) q^ + 22bqU~2 



8/5 



V = qjql 



H2 = IGpiipiqi - P292) + 4^19192 (Pi9i - 2^292) + 92 {P2 - 9i92) > 



V = qfq^ ~\ H2 = 2p\{piqi - P292) - 92 ' (l3p?9? " 80piP29i92 + ^^vlql) + 64<7?92' 

and six order polynomial integrals of motion 



2 5 

3 „ 6 



V = qi '92 

3 

V = ql q^ 



H2 = 2q^ 3 (p^^^ _ ^2^2) g| - 2piql {piqi - P292) + 9i ^92 ^ > 



3 _ 7 
2 



3/2 4 o^2^2 



5/2 



-f^2 = 9iPi - P292P1 



59i''Pi , 39^?"^ , 3 gipip2 P2 gi^ 

49I 4gi 8oi^/'' 



292^/^ 



+ 



29i 



+ 



Other examples may be found in [27J. It will be interesting to find a generic expression for 
all second integrals of motion H2 associated with the potential matrix (j3.15p . 

Case 4: Let us consider a non-degenerate geodesic matrix 



n(4) = - 



apip2 
a — b 



a-6 2 \ 



■P2 



\ 2 



■ PI hpip2 



such that T = (a + b)pip2 ■ 



(3.16) 



This matrix 11^'*) is compatible with two different potential matrices A. The first matrix is 
defined by the entry 

Ai2 = 9r/'+'<^>(929r/'). 
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For instance, if ^{z) = z'^ then 



/ ci{a + 2ad-b) -fii^^ cib{d + 1) ^ i_20i+2) 



A{4) 



2(a + 6) ^2 ^1 ^TT^^^^i 
V ' bia + b) ^2 91 2{a + b) 9i 



= ci 92^+1 . (3.17) 

It is easy to see that traceA^^ = Hf and that recursion operator produces the Hamiltonian 
only. 

In order to get the desired functionally independent integral of motion H2 we have to 
solve the equation ()3.9p . As above, there exist integrals H2 of second, forth or sixth order 
in momenta. Moreover, we sometimes have two independent solutions of the equation (j3.9p 
associated with different x's. For instance, integrable system with potentials 

T/ -2/3 -7/3 

y = qi q2 > 



has fourth order polynomial solution of (j3.9p at x = 1 

rri^) _ Pijpiqi -P2q2f _ i^phi - 44pip2gig2 + 4p^g2 , -1/3 -8/3 

2 - 2 i« 2/3 4/3 ^2 

and sixth order polynomial solution at x = 2 

^3"^^ = 4:Pi{Piqi -P2q2)^ + qi^^^q2^^^ {Wpfql -I6pip2qiq2+plq2) 

- ^qfV^'pi{pm-P2q2){2plql - Qp,P2q,q2 + plql) " A • 

91 ^2 

It can be a first example of a two-dimensional super integrable system with second, forth and 
sixth order polynomial integrals of motion, which form non-trivial Poisson algebra |27j . 

The second potential matrix A, which is compatible with geodesic matrix 11^^^ , looks like 



^ 92 1\ TT-i, — ^^ti I2 +C2gi 

2 a + b g a + b 2 o 2b 2a 
a ei „a^ „ ^2 ei(a +4a6-i) ) „a-ft „ a-6 2ac2 9i 92 
p-92 91 -^2 92 <?2 9l ^ 



2a a + b 3b— a 



A(4') 



^ - 92 9i - 9192. (3.18) 

In this case the second integral of motion may be obtained from the recursion operator 
i/^' ^ = trace _ ^2 ^ _ (g.p^ _ q.p^f _ ^^^^ ^ q^" . 
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3.3 Integrals of motion via variables of separation 

Substituting known separation coordinates u = (ui, . . . ,Vn) (|2.5p and momenta v = (vi, . . . , Vn) 
into separated relations (|1.5p and solving the resulting equations with respect to i^i, . . . , Hn 
one gets a lot of separable Hamiltonians. The main problem is to propose an effective proce- 
dure of selection natural Hamiltonians similar to the Benenti recursion procedure. 

The interim problem is to find the momenta v = (vi, . . . ,Vn) canonically conjugated 
to coordinates u = (ui, . . . ,Vn) (|2.5|) . Different algorithms for explicit computation of the 
Darboux-Nijenhuis variables {u,v) are discussed in \T9 \ [52 l [35 l H6] . 



Case 5: Let us consider geodesic matrix 

n(5) = 







(3.19) 



apip2 P2 

which may be endowed with two compatible potential matrices 



A(5) 



/ ±- 







^ ciq^ + C2(?i 







a + 2 
' a-2 



and A(^') 

In first case, separation coordinates ui^2 are the roots of polynomial 



S(5)(A) = (A - ni)(A -U2) = X^-2 (pi + cig" 



4 
a — 2 



whereas momenta vi^2 = ^(A = 141,2) are defined by the polynomial 



1 



pI - ciq, 




a + 2 
, a — 2 



q2 {a-2)q^-'pi\ q2{r + 8pi) 



2ci 



4p2 



+ (a - 2)pi 



a + 2 
. a-2 



8ci 



+ qi 



where r = trace iV. This first order polynomial A(X) is the solution of the auxiliary equations 

{^(A), Bi^,)} = , {AiX),Aif,)} = , (3.20) 

A — /X 

which ensure that values of A{X) in uj ()2.5p are the desired momenta 

Vj = A{uj) , {ui,Vj} = 6ij , {vi,Vj} = . 

Associated with matrix A*-^'-* separation coordinates ui^2 are defined by another second order 
polynomial 

B(5')(a) = (A--ui)(A-n2) = A2 

for which solution of the equations ()3.20p looks like 

/ 

^(A) ' ' 



2_ 

2P2 + Cl^i + 02^1 



4 

' a-2 



_2 

2ciC2gi +4cip^ + cl 



X 



q2 2c2 + ciq^ 



2 

a-2 



V 



2^92 



X+P2 [P2 + ciq^ 
\ 
J 



2 

a-2 



(a - 2)qr'pi 



1 / g 2(4c2 + ciqr')T + 2ci(c2gi + ci) a-2 f±| 



2 

a — 2 
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In both cases after substituting Pi^2{u,v) into common geodesic Hamiltonian T one gets 
separable Hamiltonians, which have the so-caUed Stackel form in {u, i;)-variables [6l 1191 [22t 135]. 
It allows us to easily find all the separable potentials V{ui,U2) and to construct additional 
Stackel integrals of motion. For instance, in first case, if a = —1, we have 

(5) 32c/^ Uivj + U2V^ 1 pI+pI Ci -2/3 /o 2 , a 2\ 

HI = — — = + (^1 + + U2) = — - — + -r^Qi + 492) . 

9 ^ui — Ju2 2 18 



In second case, if a = —1 and c\ = 0, we have 

o 3/2 2 3/2,9 3/2-1 2,2 
(5') %C{ V{u{ +f|U2' 1 H+P2 C2 -2/3. „ 2 , n 2^ 

^1 = + -M\ + ^^2) = — - — + ^ 9i (2^2 + 9gi) • 

9 ui — U2 4 1 6b 

If a is arbitrary, then we can obtain the same Hamiltonians after some additional canonical 
transformation of (p, q) variables. 

These bi-integrable systems are the so-called first and second Holt-like systems p!7]. The 
third known Holt-like system may be obtained from the Henon-Heiles system (|3.3p . 

In generic case the separable geodesic Hamiltonian T has more complicated non-Stackel 
form and, therefore, usually we do not know how to get separable potentials and additional 
integrals of motion. Examples of such generic Hamiltonians may be found in next Section. 



4 Bi-integrable systems on sphere 

In this Section we consider geodesic matrices on a cotangent bundle T*W^ of the sphere W^: 



I 



P' 



jk 



k=l 



dpi 



Yikiq) 



an 



ik 



dpj 



\ 



-H 



V 



E 



k=l 



dU 



ki 



an 



kj 



Zk{p) 



(4.1) 



By definition P!^ is the Poisson bivector compatible with canonical ones, so that 

[P,P^] = [P^,P^]=0. 



(4.2) 



It means that the geodesic matrix H and the functions Xjfc(g), yjfc(q), "Zk^p) are the solutions 
of these equations. As above, we restrict ourselves by particular solutions only. 

We present some examples of bi-integrable natural systems on two-dimensional sphere , 
which are related to the rigid body dynamics. In order to submit the Hamilton function into 
the standard form, we will use the vector of angular momentum J = (Ji, J2, J3) and the unit 
Poisson vector x = (xi,a;2,X3), see 

If the square integral of motion p^ = (x, J) = is equal to zero, the rigid body dynamics 
may be restricted on the sphere S^. There exists a standard spherical coordinate system on the 
cotangent bundle T*S^, which consists of Euler angles (pjO and the corresponding momenta 

9 = (91,92) = (0,6') and p = {p^^p2) = (jp^,pg) . 

defined by 



Xl 



sm sm I 



X2 = cos (p sm ( 



2:3 



cos ( 



sm cos ( 



sin( 



■p^ - cos (pPB ■ 



J2 



COS COS I 



smf 



■p^ + sm<j)p0 , 



J3 
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We use these variables in definitions of geodesic Poisson bivectors ()4.ip and potential parts A 
compatible with them. 



4.1 The Kowalevski top and Chaplygin system 

Firstly we consider a geodesic bivector Pip ()4.ip determined by the degenerate matrix H 



n(6) 

and by functions 



sm y cos^ 



a 



a G 



(4.3) 



\ cos^ Op"^ + sin^ 6'i 



yi2 = cos^^sinS + ax22(^) cos0^ , zi^2iP4>,Pe) = 0. 

Other functions Xj^ and jik are arbitrary. This matrix II^^) is consistent with the following 
potential matrix 

/ 



A(6) 



a cos ac 



b sin aip (^a sin acp — b cos acp^ cot 6 ^ 



V 



a sin a(j) — b cos ad) ] tan 



-a cos + 6sinQ<; 



a,b & 



I 



The eigenvalues of the corresponding recursion operator are the variables of separation, 
which have been studied in |45t I46|. At a = 1, 2 the separable natural Hamiltonians 



1 + 



1 



sm 



+ Pfl + 2 ( a cos a<i) — b sin a<i) ) sin' 



(4.4) 



have thoroughly familiar forms in physical variables 



(6) 



Jl + J| + 2 J| + 2ax2 - 26x1 



Uf^ = Jf + J| + 2J| + 2a(x^ - xj) - 46x1X2 . 

It is easy to see that they coincide with the Hamilton functions for the Kowalevski top and 
the Chaplygin system, respectively, see [9] and references within. 

The additional integrals of motion are fourth order polynomials in momenta: 



(4.5) 



H. 



(6) 



sm 



a cos I 



p^Pg — 2(a sin + 6 cos (p) cos P6»P(/> + sin 6 (a sin iji + 5 cos i 



b sin ( 



sm ( 



and 



(6') 



sm 



— (a cos 2(/) — 6 sin 2(/)) 



2 2 2(a sin 2(/> + 6cos 2i?ii) cos^ ^ 
^^^^ + ^^T^ ■ 

2-3cos2| 



, pe + sin^ 6 {a sin 2(/) + 6 cos 2(f>y 



sm 



+ (a^ + 6^) cos 20 . 



The corresponding separation relations are non-affine (non-Stackel) relations \45\ I46j. Note 
that these variables of separation are different from the famous Kowalevski variables. 
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Remark 8 Substituting the same variables of separation into other separation relations we 
can obtain different generalizations of bi-integrable Hamiltonians H^"^ (14. Sp . see [37 ^ 145 ^ [46]. 



4.2 The Goryachev-Chaplygin top 

The Goryachev-Chaplygin top is defined by the Hamilton function 



H 



(7) 



Jf + J| + 4J| + axi + 



4-3 cos^ e 



sm 



cos^ 



(4.6) 



which is in involution with the second integral of motion 



H. 



(7) 



2 Js ( Jl + J| + ^ 



C1X3J; 



1 ■ 



This system is a bi-integrable system with respect to natural bivector P' ()4.ip defined by 



/ ^2 , 4 - COs2 ^ ^2 



n(7) 



+ 



V 



sm 



\ 



C2 



A(7) 



and by the functions 



X22 — yi2 



2 cos^ 2 



cos sin a6 



a 



cos^ 





C2 



(4.7) 



cos^ 



Pk 



The corresponding separation variables and some another bi-integrable systems separable in 
these variables are discussed in 



Remark 9 In |48j we found some linear in momenta Poisson bivectors P for five integrable 
systems on the sphere §^ with cubic additional integrals of motion. In all these cases quadratic 
in momenta bivectors P' = PP~^P can be rewritten in natural form (j4.ip as well. 

Bivectors n('') (g^]) and H^^) (g^D have various n-dimensional counterparts on T*S" and 
variant analogs on the cotangent bundles of other Riemannian manifolds. However, in order 
to get interesting integrable systems with higher order integrals of motion we have to learn to 
construct separable natural Hamiltonians directly from the variables of separation in Stackel 
and non-Stackel cases. 



5 Generalized natural Poisson bivectors on M^" 

In this section we consider some generalizations of natural Poisson bivector (II. 9p on M^", 
which are related to various modifications of the geodesic bivector and with "potential" parts 
depending on momenta. 
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5.1 Geodesic matrix 11 depending on coordinates 

We suppose that 11 depends on coordinates q and momenta p and geodesic bivector on 
]^2n jg g^Yg]2 i^y formulae (|4.ip . For instance, let us consider degenerate matrix 



n(8) 




(5.1) 



which is the similar to matrix 11*^^^ (14. 3p associated with the Kowalevski top and Chaplygin 
system. Solving the equations (j4.2p with respect to functions ^j,kiQ),Yi,k{<l) and Zfc(p) one 
gets 

Xji = Yii = -qi , zi,2 = . 
In generic case potential matrix A is compatible with II^^) if and only if 

A22 = diqj + d2^2 +d3, dk 



(5.2) 



If di = ds 



and d2 = — ci/4, then 
/ -2cig2 - C2 



A(8) 



4cigi - 4g^ ^(3cig| + 202^2 + C3) \ 



V 



ciq2 



and the integrals of motion for the Henon-Heiles system h[^^ = T + V^^^ (j3.3p and h!^^^ (j3.5 
are in involution with respect to the corresponding Poisson bracket. 
If di = —ci and ^2 = ^3 = 0, then 



(5.3) 



A(8') 



V 



C3 - Cl 



1 + 2,1 



ciqiq2 



-ciqiq2 



2i2ciql - csql + C2) \ 



Q1Q2 



(5.4) 



and the integrals of motion for the system with quartic potential h[ = T + V^-^'^ (|3.4p and 
are in involution with respect to the corresponding Poisson bracket. 
In both cases eigenvalues ui^2 of recursion operators are the variables of separation 
for the Henon-Heiles system and the system with quartic potential, respectively. The corre- 
sponding separated relations are the standard affine Stackel relations, see |331 134j. 

Remark 10 These variables of separation ui^2 have been introduced in p3] with the help 
of the singular Painleve' expansions of the solutions of the equations of motion. Different 
properties of these variables of separation and the corresponding rational Poisson bivectors 
have been studied in [Ml EH] . 



5.2 Potential matrix A depending on momenta 

Open relativistic Toda lattice associated with An root system is an integrable system on M^" 
with following first two Hamiltonians: 

n "1 

Hi = ^^Ci + di, H2 = '^(^2^Ci + dif + Ci-i{ci + difj , 



i=l 



i=l 
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where 



Ci = exp{qi - Qi+i + Pi 



-oo, 



Jn+1 



+00 . 



This system is also known to be bi-Hamiltonian with respect to second Poisson bracket on 

Tn>2n-1 



{ck,dk}' = Ck, {ck,dk+iy = -Ck , {dk,dk+i}' = Ck 



(5.5) 



which was found in [31j. Of course, it is non-natural system in M^", but the corresponding 
Poisson bivector may be rewritten in the generalized natural form, if we introduce constant 
matrix 

/ I 1 ••• l\ /O ••• o\ 



E = a 



1 



1 

1 1 



1 







a € 



\i ... 1 0/ 



and antisymmetric matrix A with entries = e'^* so 

A = -EA, 

similar to the non-relativistic Toda lattice (j2.10p . However, let us emphasize, that in this case 
potential matrix A depends on coordinates and momenta. If we put 



n = diag(exp(-pi), . . . ,exp(-p„)) , {dpU)ij = (a - 1) ^ 



{dpA),j = (a - 1) 2^ — - + 



dpk dpj 



k=l 

dA 



dUik _ 
dpk 

kj dAki 



dU 



a- 



■33 



dpj 



dqt dqj 



(5.6) 



k=l ' k=l 

then P' is a sum of the geodesic Poisson bivector and potential Poisson bivector 



dpA A 























) 





(5.7) 



-A ' dgA 



As above, it is the Poisson bivector at A = and A 7^ 0. At any a bivector (j5.7p reduces to 
the known Poisson bivector (j5.5p on M^*^"^ and all the integrals of motion are the traces of 
powers of the recursion operator N = P{P')~^ I36j . 

Remark 11 In the similar manner as for non-relativistic Toda lattice we can get natural 
Poisson bivectors P' for relativistic Toda lattice associated with other root systems. 



5.3 Additive deformations 

Using canonical transformations and deformations of separated relations we can study more 
complicated integrable systems. For instance, let us consider a generalized Henon-Heiles 
system [n\ \^2U\ with the potential 

^ = I ,2{3q! + 16gi) + c, (2gi + I) + I + I • 
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and the second integral of motion 



Ho — H, 



(1) 



44 



C5(4p?gi + 3ciq2qf + C2qf + 8C4) ^ c^j^pjql + Ciq2qf + 4c4) 



qf 



if 



Here ffg^^ is given by (13. 5p at C3 = . 

Note, that we do not have any additional information for this system, as the Lax matrices, 
r-matrices or relations with soliton equations. Nevertheless, it is easy to directly prove that 
these integrals of motion are in bi-involution with respect to the Poisson bracket associated 
with additive deformation of the natural Poisson bivector P' (13.31) 



P = P' + 



/ 








PI y 




* 









* 


* 





V 


* 


* 


* 



-C5 



2ql 



3ci 2 






Qciql 




This additive deformation may be obtained from the natural Poisson bivector P' ()3.3p using 
trivial canonical transformation 



Pi Pi + /(gi) , where f{qi 
For generalized system with quartic potential [20] 
V 



(5.8) 



'\qt + 6qhl + 8ql) + '-^{q! + 4ql) + 



C4 C5 
1 — - -I — - 



2C3 

4 2 q.2 01 Qi 

we have to shift known natural bivector P' (|3.4p on the similar additional term 

/ 

a/— 2C5 * U P2 





p = p' + 



1! 



Pi 

* 

* * 



-C5 
2qf 



4C3 



V 



8cig| + {Sciql + 4c2)g2 5- 

92 





associated with the same trivial canonical transformation (15. Sp . 

We can apply this transformations to the Poisson bivectors P' defined by geodesic matrix 
nC^) (fSTT]) and potential matrices A^*^) ([53]) and A^^') (f5^ too. In both cases the "shifted" 
variables of separation 

ui,2 = ui^2(pi ^ Pi + fiqi)^ , (5.9) 

are defined by the initial variables ui^2 obtained in [33] at C5 = 0. It is easy to prove that the 
corresponding separated relations are non-affine in Hi^2, i-e. they are non-Stackel relations, 
similar to the Kowalevski top and the generalized Chaplygin system [45 [ 06] . These separation 
relations will be discussed in the forthcoming publication. 

Remark 12 Another example of such additive deformations of the natural Poisson bivectors 
on the plane and the corresponding separation variables may be found in 
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6 Conclusion 



We address the problem of construction of natural integrable systems on Riemannian mani- 
folds Q within the theoretical scheme of bi-Hamiltonian geometry and introduce the concept of 
natural Poisson bivectors, which generalizes the Benenti construction of the Poisson bivectors 
via conformal Killing tensors of gradient type on Q. We suppose that the proposed construc- 
tion allows us to describe a majority of known integrable systems with higher order integrals 
of motion and to find the corresponding variables of separation in common framework. 

A lot of known and new integrable systems on plane and on sphere is discussed in detail. 
These examples may be useful for creating a geometrically invariant theory, which takes the 
constructive answers to the main open questions: 

• how to get and classify all the natural Poisson bivectors P' on T*Q] 

• how to describe all the natural Hamilton functions associated with a given P' . 

Now we have some particular answers obtained by direct tedious computations only. 

We hope that the theory of natural Poisson bivectors allows us to investigate known n- 
body systems and low-dimensional exotic systems, systems with known Lax representation 
and systems without it, systems associated with the Killing tensors at 11 = and systems with 
higher order integrals of motion at 11 7^ and so on. But, of course, this theory can not become 
a universal panacea and we briefly discuss some possible generalizations and modifications of 
the natural Poisson bivectors in last Section. 

The author thanks referees for very useful suggestions. 
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